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Abstract 

One of the greatest challenges in cosmology today is to determine the nature of dark energy, 
the sourse of the observed present acceleration of the Universe. Besides the vacuum energy, 
various dark energy models have been suggested. The Friedmann - Robertson - Walker (FRW) 
spacetime plays an important role in modern cosmology. In particular, the most popular 
models of dark energy work in the FRW spacetime. In this work, a new class of integrable 
FRW cosmological models is presented. These models induced by the well-known Painleve 
equations. Some nonintegrable FRW models are also considered. These last models are 
constructed with the help of Pinney, Schrodinger and hypergeometric equations. Scalar field 
description and two-dimensional generalizations of some cosmological models are presented. 
Finally some integrable and nonintegrable F{R) and F{G) gravity models are constructed. 
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1 Introduction 

Recent measurements of redshift and luminosity-distance relations of type la supernovae indicate 
that the expansion of the universe is accelerating [I]- [2]. This appears to be in strong disagreement 
with the standard picture of a matter dominated universe. These observations can be accommo- 
dated theoretically by postulating that certain exotic matter (dark energy) with negative pressure 
domainates the present epoch of our universe. Understanding the nature of dark energy and its 
many related problems is one of the most formidable challenges in modern cosmology (see e.g. 
0-[Il])- Almost all models of dark energy work in the FRW spacetime that is in other words are 
the FRW models. Practically, all FRW cosmological models of dark energy face some difficults 
related with some nontrivial problems like the coincidence problem, the fine-turn problem and so 
on. Such and other problems of modern cosmology demand more carefully investigate the physical 
and mathematical nature of the FRW models. In this context, it is important to study integrable 
(and nonintegrable) cases of the FRW models. Here we mention the well-known fact that Ein- 
stein's gravitational equations and their some generalizations admit integrable reductions (see e.g. 
Refs. |13|-|17j and references therein). In this work, we consider some examples of integrable and 
nonintegrable FRW cosmological models induced by some known linear and nonlinear second-order 
ordinary differential equations (ODE) (see also 18 - 20 ). 

The paper is organized as follows. In section 2, we give some main informations on the FRW 
cosmological model for the flat spacetime case. In the next section, we study some nonintegrable 
FRW cosmological models. A new class of integrable FRW models was considered in section 4. 
These results were generalized for the modified f{R) and F{G) gravity models in sections 5 and 
6. Sections 7 and 8 are devoted to the scalar field description and two-dimensional generalizations 
of some models. In the last section, we give the conclution. 
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2 FRW cosmology 

We start with the standard gravitational action 

S = J V^d^x{R + Lm). 



(2.1) 



Here R is the scalar curvature and L„j is the Lagrangian of the matter. Now we consider the FRW 
spacetime with the scale factor a{t) and metric 



ds"^ = -di^ + a^{t) 



1 — fcr^ 



(2.2) 



Here k can take any value but it is related to ( — ,0,+) curvatures according to sign. One sets 
a = da/dt and computes Christoffel symbols to be 



kr 



kr"^ ' 



r^2 = adr^, = adr sin^ 6 



^ 01 



T^2 T^3 

^ 02 > ^03 



^ 22 



p2 
^ 12 



■p3 _ -1 p2 



They lead to 



with Ricci scalar 



-Roo 



33 



r(l - kr^), = -r(l - fcr^) sin^ 6», 

= — sin cos 6, = cot ^, 



(2.3) 



R 



ad + 2d2 + 2k 
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a 1 — fcr^ ' 

iJ22 = [aa + Ta^ + 2fc] , i?33 = [aa + 2a^ + 2fc] sin^ 6» 



The Friedmann equations read as 



i? = 6 



47rG 



(p + 3p), 



(2.4) 
(2.5) 

(2.6) 



where we note that p < — p/3 implies repulsive gravitation. Recaling the Hubble parameter 
B. = da~^ these equations can be written as 



H' = ^p-^, H = -4nG{p + 3p) + ^. 

fj Qi CI 



(2.7) 



The Friedman equations with a cosmological constant for a homogeneous isotropic universe models 
are 



A k 8ttG d A 47rG, 
2+^~^' 7 = T —{P + ^P), 



a 



(2.8) 



In this work, wc consider the case: fc = A = and set SttG =1. If the FRW spacetime is filled 
with a fluid of energy density p and pressure p, then conservation of energy-momentum tensor 



VT^. = 0, 



(2.9) 



gives 



p + 3Hip + p) = 0, (2.10) 

where a dot represents differentiation with respect to t. So finally the equations for the action 
(2.1) we can write in the H-form 

t2 



p = -2H - 3H^ 
p = -3H{p + p) 



(2.11) 
(2.12) 
(2.13) 
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or in the N-form 



p = -27V -3iV^ (2.14) 
p = 3iV^ (2.15) 
p = -3N{p + p), (2.16) 

where N = In a. In this work, we study some FRW cosmogical models, for example, with the 
equation of state (EoS) [3]-[2 (see e.g. also [6]-[T2]) 

V = /i(A^, ^, t)p + /2(7V, TV, t)/ + /3(iV, TV, t), (2.17) 

where /3 7^ 1 and /i = /i(TV, iV, i) are some functions of TV, iV and t. 

3 Nonintegrable FRW cosmological models 

In this section, we consider some known and new FRW models, for which, iV satisfies some second- 
order hnear and nonhnear ODEs. These ODEs are nonintegrable so that the corresponding FRW 
cosmological models are nonintegrable (see e.g. [TSJ-IIS])- Before we go any further, we explain 
more clearly (as mathematical claim) what is integrable FRW cosmology. If one of dependent 
variables of some FRW model, say N for our case, satisfies some integrable differential equation then 
such FRW model we call an integrable FRW model (or in other words integrable FRW cosmology) . 
On the other hand, if N satisfies some nonintegrable differential equation then such FRW model 
we call a nonintegrable FRW model (or in other words nonintegrable FRW cosmology) . Now let us 
give the definition of integrable equations. In mathematics and physics, there are various distinct 
notions that are referred to under the name of integrable systems (equations) . In our case, the given 
nonlinear differential equation is said to be integrable if it admits infinite number of independent 
conserved quantities (that is integrals of motion) (/i, /2, ...) in involution. Two conserved quantities 
are said to be in involution if they Poisson-commute that is {Ii,Ij} = where {/,<?} is 
some Poisson bracket. 

3.1 ACDM cosmology 

We start with the ACDM cosmology. Let the parametric EoS has the form 

p==-A, p^3N^. (3.1) 

It is well-known that in this case, N satisfies the following equation 

TV = 0.5A- 1.5TV2. (3.2) 

The corresponding EoS parameter reads as 

uj = -l + ^V^. (3.3) 

3.2 Pinney cosmology 

In this subsection, we study the cosmological models induced by the Pinney equation. This equation 
we here write as [22] - [23] (see also e.g. [39] ) 



v = my + ^, (3.4) 



yi 



where ^ — ^(t), k — const. 

1) Let the parametric EoS has the form 



p ^ -iN'^ + 2(^{t)N ~2kN-^, (3.5) 
p = 3TV2. (3.6) 
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Then it can be shown that N satisfies the Pinney equation 

N = -em+^, (3.7) 

where C = C(^)i ^ — const. It is known that as C = 1, the Pinney equation has the following 
particular solution (see e.g. |21| ) 



Then 



N = (cos^ t + fcsin^ i) • (3-8) 



N = {k-l)smtcost{cos'^t + ksm'^ty°-^, (3.9) 
N = [A:- (cos'*i + 2fcsin^icos^i + fc^sin*i)](cos^t + A:sin^t)"^-^ (3.10) 



and 



3(fc - l)2sin^tcos2i 
cos^ t + k sin^ t 

3(A: - 1)2 sin^ t cos^ t 2[k - (cos^ t + 2k sin^ t cos^ t + P sin" t)] 



= , 2^ ' (3-11) 



P = ^ 2 ^ 2 ■ (3-12) 

cos^t + Zcsin t (cos^i + fcsin ty-^ 

The last equation is the parametric EoS corresponding to the Pinney equation (3.7) [if exactly to 
its solution (3.8)] which we can write in the usual form as 

p= -p-2{kN'^ - N). (3.13) 

The corresponding EoS parameter takes the form 

1 2k 
3 ~ 3(cos2t + fcsin^ ty 



2) Let us we consider another derivation of the Pinney cosmology. For the flat FRW metric 
(2.2), the usual Einstein-Dirac equation has the form 

3H^-p = 0, (3.15) 
2H + 3H^+p = 0, (3.16) 
ij; + 1.5Hij + tj°V^ = 0, (3.17) 

iP + l.bHijj - iV^j° = 0, (3.18) 
p + 3H{p + p) = 0, (3.19) 

where the kinetic term, the energy density and the pressure take the form 

y = 0.5i(Vi7V - ^7V), P^V, p = uVu-V, u = i)^. (3.20) 
Let the potential has the form 

V ^Vo + mu + ici^iu^i - 3c~^^2W^ - l.Sc^^^au^, (3.21) 
where — S,j{t) are some function of t and Vb, m,c = consts. In this case we have 

Y + 3Hu{Vu + uVuu) = 0, V + 3HY = 0, u = ca-^. (3.22) 

As 



Y 
2a 

for the potential (3.21), we get that the scale factor satisfies the following equation 

ann = Uv)a + 6(^)a"' + Uivh'^ + Uv)a-\ (3.24) 
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where drj = adt, ^4 = 0.5cm. Consider particular cases. 

i) Let ^1 = —d'^{ri), ^2 ~ k, ^4 = ^5 = 0. Then the equation (3.24) becomes 

ar,n + 0^iv)a~ ^ ^0. (3.25) 

It is the Pinney or Ermakov-Pinney equation. The general solution of Eq. (3.25) is given by [2 3) 



a= y^Az2 + Bz| + 2CziZ2, (3.26) 
where Zj (77) are linearly independent solutions to the equation 

zj^^ + 9^{rj)zj ^0. (3.27) 
Here the constants A, B, C satisfy the constraint 

AB-C'^ = kW-'^ (3.28) 

and 

W = ZiZ2jj - Z2Zxr, (3.29) 

is the Wronskian. 

ii) We now consider the case: = — = const, £,2 — £,4 — 0, ^5 = — k. Then Eq. (3.24) 
takes the form [33] 

arin + 9la+^ = 0. (3.30) 

3.3 Schrodinger cosmology 

We now reconstruct the FRW cosmological model induced by the linear Schrodinger equation 

N^uN + kN, (3.31) 

where u = u{t), k ~ const. Then the parametric EoS takes the form 

p = -3N^ -2uN -2kN, (3.32) 

p = 3N^. (3.33) 

It is well-known that as m = n{n — l)t^^ — k, Eq.(3.31) has the following particular solution 

TV = At". (3.34) 

Then 

p = 3A2n2t2("-i), (3.35) 

p = -3A2n2i2(n~i)„2n(n-l)At"-2 (3.36) 

or 

p= -p-2n{n~l)X'^3^^^p^^^. (3.37) 
The corresponding EoS parameter takes the form 
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3.4 Hypergeometric cosmology 

i) We start with the equation 

N = t-^{l - t)~^{[{a + b+l)t- c]N + abN}. (3.39) 

It is the hypergeometric differential equation. In this case, the parametric EoS has the form 

p = -3N'^-2t-^{l-t)-^{[{a + b+l)t-c]N + abN}, (3.40) 
p = 3Ar2. (^34^-) 

The solution of the equation (3.39) is the hypergeometric function 

N = 2Fi{a,b;c;t). (3.42) 

Some particular solutions are: 

ln(l + t) = t2Fi(l,l;2;-t), (3.43) 

{l-ty = 2F,{a,b;b;t), (3.44) 

arcsint = t2i^i(0.5,0.5; 1.5;*^), (3.45) 

and so on. As an example, let us consider the solution (3.45) that is N = arcsinf. Then 

N = {l-t^)-°-^, iV = t(l -f2)-i-5. (3.46) 



and 



Hence we obtain 



or 



p = 3il-t^)-\ (3.47) 
p = -3(l-t^)-^ -2t(l-i2)-i-5. (3.48) 

p= -p-3-^-^2tp^-^ (3.49) 



p = -p-3-^-^2{p-3)°-^p. (3.50) 
For this case, the EoS parameter takes the form 

oj = -l- ^^=. (3.51) 

ii) As the next example of the hypergeometric cosmology we can consider models induced by 
elliptic integrals. As an example, let us consider the complete elliptic integral of the first kind K: 

Assuming A'' = K{t), we can calculate all expressions for this case. But we drop it as this case is 
the particular reduction of the model (3.42) due to of iV = K{t) = 0.57r 2-Fi(0.5, 0.5; 1; t). 

iii) The last example is the case when N is equal to one of incomplete elliptic integrals. For 
example, we can put 



m-Fm^J^ ^^_,^_^^^^ ^ (3,53, 



Similarly to the previous cases we can also find all expressions to describe the cosmological model 
but omit it. 



7 



3.5 Weierstrass cosmology 
3.5.1 The Weierstrass gas 

One of the most interesting examples of Weierstrass cosmologies is the Weierstrass gas which has 
the following EoS 

p=-B[p[p)f-'>. (3.54) 



The well-known Chaplygin gas model 



B 

p= (3.55 

P 



is the limit (or particular case) of the Weierstrass gas when the Weierstrass function takes the form 

p{p)^p-^ (3.56) 

3.5.2 The generalized Weierstrass gas 

The generalized Weierstrass gas corresponds to the EoS 

p = -i3[p(p)]°-^". (3.57) 
Its limit (or the particular case) is the generalized Chaplygin gas [39] 

P=-^ (3.58) 

that follows from (3.56) and (3.57). 

3.5.3 The modified Weierstrass gas 

At last we present the modified Weierstrass gas (MWG) model. Its EoS is given by 

p-Ap-i3[p(p)]0-5". (3.59) 
The MWG is some generalization of the modified Chaplygin gas [ID] 

p = Ap-^. (3.60) 

Finally we would like to give the more general form of the MWG. Its EoS reads as 

p = A[p(p)]-°-5-S[p(p)]"-^". (3.61) 

4 Integrable FRW cosmological models 

After the previous section, we think that logically we are here in the position to make the next 
step, namely, to consider some integrable FRW cosmologies (see e.g. [IH]-[in])- Here we restrict 
ourselves to Painleve cosmology that is cosmologies induced by the well-known Painleve equations. 
These six Painleve equations (P/ - Pv/ equations) were first discovered about a hundred years ago 
by Painleve and his colleagues in an investigation of nonlinear second-order ordinary differential 
equations. Recently, there has been considerable interest in the Painleve equations primarily due 
to the fact that they arise as reductions of the soliton equations which are solvable by Inverse Scat- 
tering Method (see e.g. [26] V Consequently, the Painleve equations can be regarded as completely 
integrable equations and possess solutions which can be expressed in terms of solutions of linear 
integral equations, despite being nonlinear equations. The Painleve equations may be thought of a 
nonlinear analogues of the classical special functions. They possess hierarchies of rational solutions 
and one-parameter families of solutions expressible in terms of the classical special functions, for 
special values of the parameters. Further the Painleve equations admit symmetries under affine 
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Weyl groups which are related to the associated Backlund transformations. Although first discov- 
ered from strictly mathematical considerations, the Painleve equations have arisen in a variety of 
important physical applications including statistical mechanics, plasma physics, nonlinear waves, 
quantum gravity, quantum field theory, general relativity, nonlinear optics and fibre optics. In this 
section, we show that these six Painleve equations also can describe, in particular, the accelerated 
expansion of the universe that is the dark energy. 

4.1 Painleve cosmology 

We work with the N-forni of the Einstein equations that is with the system (2.14)-(2.I6). We now 
assume that N satisfies one of Painleve equations. Consider examples [below a, /3, 7, d, k and /i are 
arbitrary constants]. 

4.1.1 P/ - model 

Let the parametric EoS has the form 

p = -3N^ _ 12N'^ ~ 2t, (4.1) 
p = 3N^. (4.2) 

Then we can show that N satisfies the P/ - equation [25] 

N = 6N^+t. (4.3) 

It is well-known that the P/ - equation can be expressed as the compatibility condition 

At -Bx + [A, B]=0 (4.4) 

of the following linear system 

ox - (4-^) 
9$ 

where A = A{t, A) and B = B{t, A) are matrices, A is a spectral parameter which is independent 
of t. For the P/ - equation, A and B have the form 

A = (4A'^ + 27V2+t)CT3-i(4A^A^ + 2iV2 + i)cr2-(2AiV + 0.5A"i)o-i, (4.7) 
B = {\ + \-^N)as-i\-^Na2, (4.8) 

where Pauli matrices are given by 

l\ /O -A fl 



4.1.2 P// - model 

Our next example is the case when N satisfies the P77 - equation 

N = 2N^ + {t- to)N + a, (4.9) 

where N = N(t, a), a = const. Then the parametric EoS has the form 

p = -^m^ -4N^ -2{t~to)N -2a, (4.10) 
p = 3N^. (4.11) 

We note that the P// - equation (4.9) can written as the compatibility condition (4.4) of the linear 
system (4.5) - (4.6) with 

A = -i{4:\^ + 2N^ + t)<73 - 2N<72 + i4:\N - a\''^)ai, (4.12) 
B = -iXa3 + Nai. (4.13) 
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It is well-known that the P// - equation has the following particular solution 

N = Nit,l)^-it-to)-'. 

Then 



-3 



and 



P = 3{t- ^o)"^ P = ^t- to)-' - Ht - to)-\ 

For this example, the EoS and its parameter take the form 

, 0.75 



P 



and 



-1+3(^-^0). 



(4.14) 

(4.15) 
(4.16) 

(4.17) 
(4.18) 



So if t < io + 0.5 {t > to + 0.5) then this P//-model describes the accelerated (decelerated) 
expansion of the universe and the case t = to corresponds to the cosmological constant. 

4.1.3 Pin - model 

Let the parametric EoS has the form 



p = -3iV^ - 2 
p = 5N^. 
Then N satisfies the P/// - equation [26] 



N = —N^ --(N- aN"^ - /3) + iN^ + — , 



(4.19) 
(4.20) 

(4.21) 



where N = N{t, a, /3, 7, 5). One of particular solution of the P/// - equation has the form [55] 

(4.22) 

Then 



N = N{t, a, 0, 0, -aK^) = Kt^'^ 



-r 



5/3 



and 



--4/3 , 4« .-5/3 



The corresponding EoS and its parameter take the form 

4 



P 



and 



-1 



^27^6 
4 



1.25 



(4.23) 
(4.24) 

(4.25) 
(4.26) 



4.1.4 Pjv - model 

Let us consider the parametric EoS of the form 



P 
P 



'3N^ - 2 



3iV^ 



5_ 

N 



(4.27) 
(4.28) 
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Then for N we come to the equation 

1 



N = —N^ + l.5N^ + AtN^ + 2{f - a)N + — , 
2N ' N 



(4.29) 



which nothing but the P/y - equation. Here N — N(t,a,6). One of particular sohition of this 
equation has the form |26] 

iV = 7V(t,0,-2) = -2t. 



Then 
and 



N = 



p=12, p = -l2. 
The corresponding EoS and its parameter take the form 

p = -p, u = -1. 

4.1.5 Py - model 

Let N is the solution of the Py - equation [26] 



(4.30) 
(4.31) 
(4.32) 

(4.33) 



_i, 5N{N+l) 



where N = N{t, a, /3, 7, S). The corresponding parametric EoS is given by 
-3N^ - 2 



P 
P 



37V2, 



(4.34) 

,(4.35) 
(4.36) 



where / = SN{N + 1){N — 1) ^. It is well-known that Eq.(4.34) has the following particular 
solution (see e.g. [22] and references therein) 



Then 
and 



N = N{t; 0, 0, p, -0.5p^) = ne''*. 



The corresponding EoS and its parameter take the form 



and 



4.1.6 Py/ - model 



6K 



(4.37) 

(4.38) 
(4.39) 

(4.40) 
(4.41) 



Our last example of integrable FRW cosmologies is the Py/ - model. Its the parametric EoS is 
given by 



P 
P 



-'iN^ - 2 



3iV^ 



0.5 — 



1 



1 



1 



N N -\ N-t 



1 1 



1 



t t-1 N-t 



N + J 



(4.42) 
(4.43) 
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where J = t''^{t - 1)-^N{N - 1){N - t) [a + /3tN-^ + j{t - 1){N - l)"^ ^ ^^^^ _ i^^j^ _ i)-2j 
The corresponding nonhnear ODS is the Pvi - equation 

iV = 0.5 ( 4 + ^J— + ^J—\ N'^-(- + + -J— I N 



N N-1 N-tJ \t t-l N-t 

+ t-^it - 1)-^N{N - 1){N ~t)[a + ptN~^ + j{t -1){N - 1)-^ + 5t{t - 1){N - t)'^] , (4.44) 
where N = N{t, a, l3, 7, 6). Its particular solution is [26 

N = N{t; 0.5^2, -O.Sk^, 0.5M^ 0.5(1 - /i^)) = i" ^ (4.45) 

Then 

N = 0.5^"°■^ N = -0.25*-^-^ (4.46) 

and 

P=lt-\ p=-^t-'+OM-'-'. (4.47) 
The corresponding EoS and its parameter take the form 

P = -p+^^p'-' (4.48) 



-l + lt-°-^. (4.49) 



and 



4.2 Hamiltonian structure 

It is very important that all Painleve equations can be represented as Hamiltonian systems that 
is as (see e.g. and references therein) 

dF 

q = V, (4.50) 
or 

dF 

r ^ (451) 

where F(q, r, t) is the (non-autonomous) Hamiltonian function. Consider some examples (see e.g. 
[26] and references therein). 

1) P/-model. In this case, q,r,F read as 

q = r, (4.52) 
r = Qq^ + t, (4.53) 
F = 0.5r2 - 2q^ - tq. (4.54) 

2) P//-model. In this case we have 

q = r-q^-O.bt, (4.55) 
r = 2gr + a + G.5, (4.56) 
F = 0.5r2-(g2 + o.5t)r-(a + 0.5)g. (4.57) 

3) P///-model. In this case we have 

tq = 2q^r - k2tq^ - {29i + l)q + kit, (4.58) 
tr = -2qr^ + 2k2tqr + {20i + l)r-k2{{9i+92)t, (4.59) 
tF = q^r^ -[k2tq'^ + {20i + l)q-kit]r + k2{0i+02)tq. (4.60) 



5 F{R) - gravity models induced by second-order ODEs 

In this section, we consider F{R) - gravity models induced by second-order ODSs. Some of these 
F(R) - gravity models are integrable and others are nonintegrable. 
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5.1 Integrable F{R) - gravity models 

Our aim in this subsection is to present a class of integrable F{R) - gravity models for the FRW 
metric case. To do it, we use again Painleve equations. Consider the action of F{R) - gravity [37] 
(see also e.g. [28]-[34]') 

S = J ^d*x[R + f{R) + L^]. (5.1) 
In the case of the FRW metric, the equations for the action (5.1) are given by 

3H^ = Peff, 2H + 3H^ ^-peff. (5.2) 

Here 

Peff = -0.5f + 3{H^ + H)f'-18H{H + iHH)f", (5.3) 
p^ff = 0.5f^{3H^ + H)f'+6{H + 6HH + 4H^ + 8H^H)f''++36{H + AHHff'''i5A) 

where / = df /dR etc. To construct integrable f{R) - gravity models, we assume that the function 
f{R) satisfies some integrable ODE. As an example, we here demand that f{R) is a solution of 
one of Painleve equations. Let us present these equations. 
1) Fi{R) - models, a) Fia{R) - model: 

f"=6f + R (5.5) 

or 

f" = 6F^ + R. (5.6) 
b) Fj B (R) - model. Note that instead of these two models we can consider the following ones 

/ - 6/2 + t 



or 



F^6F^ + 1. 



2) Fii{R) - models, a) Fiia{R) - model: 

f"=2f + Rf + a (5.7) 



or 

f" = 2F^ + RF + a. (5.8) 
b) Fiib{R) - model. The alternative models are given by 

/ = 2/3 + tf + a 

or 

F ^2F^ + tF + a. 
3) Fiii{R) - models, a) Fiiia{R) - model: 

/" = jf''-Jl (/' - -l3)+lf + j (5.9) 

or 

f" - -^{F - aF' - /3) + + i. (5.10) 

b) Fjjib{R) - model. The alternative models are given by 

f = jf'-]if-o.f-P)+lf + j 
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or 

4) Fiv{R) - models, a) Fiva{R) - model: 

/" = ^/'' + l-5f + 4i?/2 + 2{R' -a)f + j (5.11) 

or 

1 A 

f" = — + l.SF^ + 4RF'^ + 2(i?2 - a)F + - . (5.12) 



b) Fivb{R) - model: 



/ = + 1.5/3 + 4^/^ + 2(^2- a)/ + I 



or 

1 (5 
F = —F^ + 1.5F^ + 4tF^ + 2{f - a)F + -. 

5) Fv{R) - models, a) Fva{R) - model: 



or 



b) Fvb{R) - model: 

/ = + j^)/^ - ^(/ - 7/) + i-^(/ - Ifiaf + ^) + ^^^i:^ 

or 

F=i^ + J^)F' -^iF-lF) + r\F - l)\aF + PF'^) + 
6) Fvi{R) - models, a) Fvia{R) - model: 

1 1 1 \ ,'2 / 1 1 1 



+ R-''{R- 1)-VW- [a + + 7(fl - - 1)"" + SR{R (5-15) 

Instead of this equation we can consider the following one 



F F-l F-RJ \R R-1 F-R 

+ R-^{R - ly^FiF - 1){F -R)[a + pRF'"^ + 7(i? - 1){F - l)-^ + 5R{R - 1)(F - i?)-^] . 

(5.16) 

b) Fvib{R) - model: 



/ /-I f-tj' \t t-i f-t^ 

+t-'{t - 1)-V(/ - 1)(/ -t)[a + /3tr' + l{t - !)(/ - 1)-' + 5t{t - 1)(/ - t)-^] 
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+t-\t - l)-^F{F --1){F~ t) [a + ptF-^ + 7(t - 1)(F - l)-^ + 5t{t - 1)(F - i)-^] . 
In the above, a,/3,7, ^ are some constants and the function F{R) corresponds to the action 



S = j ^d''x[F{R) + Lra]. 



(5.17) 



Some comments in order. All above presented cosmological models [J — 

/, /J, ///, IV, V, VI) are integrable due to of integrability of Painleve equations. In particular, 
this means that these models admit n — soliton solutions. As an example, let us present here some 
exact solutions for the Fu - cosmology given by the equation (5.7). It is well-known that Eq. (5.7) 
has the following rational solution 



fiR) = f{R;n) 



d_ A f En-i{R) 
dRy^'y En{R) 



(5.18) 



where the En{R) are monic polynomials (coefficient of highest power of R is 1) and satisfy the 
following equation 



En+i{R)EN-i{R) = REliR) + ^E^{R) - AEN{R)E'^iR). 
This equation has the following first solutions ^26j 



(5.19) 



(5.20) 
(5.21) 
(5.22) 
(5.23) 
(5.24) 
(5.25) 
(5.26) 

where L ^ 1448832000i?/^ + 1931776000007?^ - 38635520000 and so on. The corresponding expres- 
sions for the function f{R) have the form ^6} 



Eo{R) 


= 1, 




Ei{R) 


= R, 




E2{R) 


^ i?^ 4 


-4, 


EsiR) 




- 20R^ - 80, 


Ei(R) 




+ 60i?^ + 11200i?, 


E^iR) 




-f 140i?^2 _^ 2800i?^ -1- 78400i?^ - 313600i?3 - 6272000, 


EeiR) 




+ 280i?^^ + 18480i?^^ + 627200i?^^ - 17248000i?^ + L, 



f{R) 


= f{R-A) 


fiR) 


^f{R-:2) 


fiR) 


= /(i?;3) 


fiR) 


= f{R-A) 



1 

^R' 

1 3i?2 



R i?3+4' 
3i?2 6i?2(i?3 _^ 10) 



i?3+4 i?6 + 20i?3 - 80' 
1 6R'^{R^ + 10) 9R^R^+40) 



R R^ + 20i?3 - 80 -f 60i?6 -I- 11200 



(5.27) 
(5.28) 
(5.29) 
(5.30) 



and so on. It is also interesting to note that the function f{R) can be expressed by the so-called 
T-function as [25] 

'T„_l(i?)^ 



where 



TniR) 



= f{R;n) 




Pi{R) 


PsiR) 


PiiR) 


PsiR) 


P^'-'\R) 


pt'\R) 



Tn{R) 



P2n-l{R) 
P2n-l{R) 



pt'-hR) 



(5.31) 



(5.32) 



Here Pj{R) are the polynomials defined by Pj{R) — for j < 0, and 



^J2p,{R)X^. 



(5.33) 
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5.2 Nonintegrable f{R) - gravity models 

Obviously that some F{R) - gravity models can be constructed by nonintegrable second-order 
ODEs. Let us consider examples. 

i) Our first example is the hypergeometric differential equation (see e.g. |28) ) 

i?(l - R)f" + [c-{a + b+ l)R]f' - abf ^ 0. (5.34) 

It has the solution f{R) = 2F{a, b; c; R) which is the hypergeometric function. We can also consider 
the i- version of the equation 

t{l - t)f + [c - (a + 6 + l)t]f -abf = (5.35) 

with the solution f{t) — 2F{a,b; c;t). 

ii) Another example is the case when ,f{R) satisfies the Pinney equation 



/"+a(i?)/ + ^=0. (5.36) 



If ^1 = 1. ^2 = = const, this equation has the following solution [5T] 

f{R) = cos^ R + K^si-a^ R. (5.37) 
Its "i-form" is /(<) = cos^ t + sin^ t which is the solution of the Pinney equation 

/ + ei(t)/+^=0 (5.38) 

as ^1 = 1, ^2 = K = const. 

iii) Let us we present one more example. Let the function / satisfies the equation 

/" = 6/2 - 0.552, (52 = const) (5.39) 

or 

/ = 6/2 - 0.552. (5.40) 
These equations admit the following solutions 

f{R) = p{R) (5.41) 

and 

fiR) EE fit) = Pit), (5.42) 
where p(i?) and pit) are the Weierstrass elliptic functions. 

6 F{G) - gravity models induced by second-order ODEs 

The next important modified gravity theory is FiG) gravity. Let us now we extend results of 
the previous section to the F{G) gravity case that is consider F(G') - gravity models induced by 
second-order ODSs. As in the previous FiR) gravity case, some of these FiG) - gravity models 
are integrable and others are nonintegrable. Consider examples. 

6.1 Integrable F^G) - gravity models 

In this subsection, we present a class of integrable FiG) - gravity models for the FRW metric case 
using again Painleve equations. The action of FiG) - gravity we write as [3T]-[32] (see also e.g. 

m-m) 

S = J V^d^x[R + /(G) + L„] (6.1) 
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or 

S = J ^/^d^x[F{G) + Lm]. (6.2) 

Here 

G = R''- AR^.R^'' + R^^^rR"""^ (6.3) 
is the Gauss-Bonnet invariant which for the FRW metric takes the form 

G = 24H'^{H + H^). (6.4) 

The FRW-equations for the action (6.1) are given by 

3H^ = PG+Pm, 2H + 3H^ = -{pG+Pm). (6.5) 

Here 

PG = GfG-f-24H^G)fGG, (6.6) 

PG = -pG + 8H^G^fGGG + E. (6.7) 

where £; = -192fGG{'iH^H -SH^HH -6H^H^ - H*H -3H^H -ISH^H^) and /g = df/dG etc. 
To construct integrable /(G) - gravity models, we assume that the function /(G) satisfies some 
integrable ODE, namely, one of Painleve equations. We here just list such integrable models, the 
investigation of which leave for the future studies. Let us present these equations. 
1) Fi{G) - models. 
a) FiAiG) - model: 

fGG = 6f + G (6.8) 

Fgg = 6F2 + G. (6.9) 
f = 6f + t 
F = 6F'^ + 1. 

fGG = 2f + Gf + a (6.10) 

Fgg = 2F^ + GF + a. (6.11) 
f = 2f + tf + a 
F = 2F^ + tF + a. 



or 



b) Fib{R) - model: 



or 



2) Fii{G) - models. 
a) Fiia{G) - model: 



or 



b) Fiib{G) - model: 



or 



3) Fiii{G) - models. 
a) Fiiia{G) - model: 



or 



b) Fiiib{G) - model: 



fGG = -^fh -^ifc- -P)+7f + j (6.12) 
Fgg = - ^{Fg - aF^ -^)+^F^ + l.. (6.13) 
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or 

4) Fiv{G) - models. 

a) Fiva{G) - model: 

fcG = Yffc + 1-5/' + + 2(G2 -a)f + j (6.14) 

or 

Fgg = ^F^ + 1.5F3 + + 2{G^ - a)F + |. (6.15) 

b) Fivb{G) - model: 

/ = + 1-5/'' + 4t/2 + 2(t2 - a)/ + I 

or 

1 A 
= -=F'^ + l.SF^ + 4iF2 + 2(i2 - a)i^ + - . 

5) Fv{G) - models. 

a) Fva{G) - model: 

fGG = + -j^)fG - ^(/g - 7/) + G-2(/ - l)\cf + + ^^yzT^ (6.16) 

or 

= (^ + p^)Fh - ^{Fg - iF) + G-\F - ifiaF + PF-') + (6.17) 

b) Fvb{G) - model: 

or 

F=i^ + J^)F' -^{F-lF) + t-\F - l)\aF + /JF"^) + 

6) Fvi{G) - models. 
a) Fvia{G) - model: 



J /-I f-GJ-"' \G G-l f-G, 
+ G-\G- 1)"2;(/- 1)(/- G) [a + PGf-^ + 7(G - 1)(/ - 1)-^ + SGiG - 1)(/ - G)-^] (6.18) 
or 



¥^f-i^f-g) [g^ g-1^ f-g^ 

+ G-'^{G - ly^FiF - 1){F -G)[a + pGF'^ + 7(G - 1)(F - l)-^ + ,5G(G - 1)(F - G^^] . 

(6.19) 

b) Fvib{G) - model: 

^t-\t - ir'fif - 1)(/ -t)[a+ pt,r^ + j(t - 1)(/ - 1)-2 + st{t - 1)(/ - 1)-'] . 

and 

V-F F-1 i^-ty \t t-i F-tJ 

+t-^{t - l)-^F{F - 1){F -t)[a + ptF^^ + 7(t - 1){F - 1)"^ + 6t(t - 1)(F - t)'^] . 

So all these models are integrable that is they admit n-soliton solutions, infinite number commuting 
integrals, Lax representations and so on. 
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6.2 Nonintegrable f{G) - gravity models 

We now briefly list some nonintegrable F{G) - gravity models which can be constructed by nonin- 
tegrable second-order ODEs. Let us consider examples. 

1) The hypergeometric differential equation (see e.g. [28] ) 

G(l - G)fGG + [c - (a + 6 + 1)G]/g - abf = (6.20) 

or 

t{l - t)f +[c-{a + b+ !)<]/ - abf = 0. (6.21) 

Note that these equations have the following solutions /(G) — 2F{a, b; c; G) and /(G) = f{t) — 
2F{a, b; c; t), respectively. 

2) The Pinney equation 

/gg + 6(G)/ + ^ = (6.22) 



or 

I 



f + + (6.23) 



If ^1 = 1, £,2 = K = const, these equations have the following solutions [21] 

/(G) = cos^ G + K^ sin^ G (6.24) 

and 



respectively. 

3) Our last models are given by 



or 

/■2 



/(G) = /(i) ^cos^t + K^sin^i, (6.25) 

fcG = 6/2 - 0.552, (52 = const) (6.26) 
/ = 6/2 - 0.552 (6.27) 



which admit the following solutions 

/(G) - p(G) (6.28) 

and 

/(G) = fit) = Pit) (6.29) 
respectively. Here p(G), pit) are the Weierstrass elliptic functions. 



7 Scalar field description 

As is well-known scalar fields play an essential role in modern cosmology since they are possible 
candidates for the role of the inflaton field driving inflation in the early universe and of the dark 
energy substance responsible for the present cosmic acceleration. A way to describe the above 
presented cosmological models from field theoretical point of view is to introduce a scalar field 4> 
and self-interacting potential U (0) with the following Lagrangian: 

= 0.502 -i7((/.). (7^^) 

The corresponding energy density and pressure are given by 

p = 0.502 + [/((/.), p = 0.502 - C/(0). (-7 2) 

Consider the P//-cosmology (subsubsection 4.1.2). Then we get that 

= / ^/ -2i2m +tN + a)dt, U = 3N^ + 2N^ + tN + a. (7.3) 
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Hence for the particular solution (4.7) we obtain 

= </,o-4t-°-5, f/ = 3r*-2r^ (7.4) 

It gives 

U = 2-163(0 _ 0^)8 _ 2-"(</) - cpof. (7.5) 

Similarly we can find the expressions of and its self-interaction potentials for the other examples 
considered above. 

8 Two-dimensional generalizations 

The above considered FRW models are one-dimensional. Here we present their two-dimensional 
generalizations following of |26j. 

a) As an example, let us consider the P//-model. We introduce a new function v{y,fi) as 

v{y,fi) = ^N{t), (8.1) 

where N{t) is the solution of the P77 - equation (4.9), y = t^fijl is a new "time" -coordinate, /i is 
some physical parameter. Here as example, we take /i = A, where A is the cosmological constant. 
Then the function v ~ v(y,A) obeys the following equation [26j 

VA — Qv'^Vy + Vyyy — 0. (8.2) 

It is the famous modified Kortewegde Vries (mKdV) equation which is integrable. 

b) Similarly, we can construct two-dimensional generalizations of the other models. For exam- 
ple, let V has the form 

v{y,A)^-^fiR), (8.3) 



/3A 

where /(i?) is the solution of the Fji - equation (5.7), y — i?v'^A. Then the function v = v{y,A) 
again satisfies the equation (8.2). Next, if we take v in the form 

v{y,A) = ~^[N'it) + N\t)], (8.4) 

V 9A^ 

where N{t) is again the solution of the P// - equation (4.9), then the function v satisfies the 
Kortewegde Vries (KdV) equation [23] 

+ 6VVy + Vyyy — 0, (8.5) 

which is also integrable. For the Fn - model (5.7), v is given by 

y{y,A) = -^[f'{R)+f{R)] (8.6) 

V 9A^ 

and y = R-^'SA, where f{R) is the solution of the equation (5.7). 

c) Now we consider the two-dimensional generalization of the equations (4.21) and (5.9). Let 
y — tfi~^ and N{t) is the solution of the P///-equation (4.21) with a = —f3 = 0.5, 7 = (5 = 0. 
Then the function 

viy,fi)^tlnNit) (8.7) 
satisfies the sine-Gordon equation [23] (see also [55] ) 



If J/ = Rfi 1 and v{y,i.i) ~ ilnf{R) then f{R) satisfies the equation (5.9). 
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9 Conclusion 



In this work we have considered some FRW cosmological models with the parametric EoS, where 
the e-folding N plays the role of the parameter. To construct the cosmological model in the explicit 
form we demand that N is the solution of some linear or nonlinear second-order ODSs. If such 
ODS are integrable then the corresponding FRW models are also integrable. It means that such 
models admit all ingredients of integrable systems such as n-soliton solutions, commuting integrals, 
Lax representations etc. 

Here we have constructed integrable FRW models induced by Painleve equations. Then we 
have extended our results for modified F{R) and F{G) gravity theories. We have obtained the 
explicit forms of f{R) and /(G) functions for a fiat FRW universe filled by some exotic fluid 
with the parametric EoS. We have then described one of considered models as FRW cosmological 
model having a scalar field and found its self-interacting potential. Note that all FRW models 
considered above are one-dimensional. Also we have discussed the two-dimensional extensions of 
some models. In this second coordinate we take some physical parameter /i, for example, 

the cosmological constant A. 

The main idea of our work is the reconstruction some gravitational models using the solutions 
of some differential equations (in our case, second-order ODSs). Here we would like to note that 
this approach works not only for second-order ODSs and can be extended to other ODSs and even 
to partial differential equations. Let us briefly demonstrate this possibility. As an example, we can 
consider the case when iV is the solution of the following equation 

^ ^ g^N ^ ga, (9.1) 

where (72, 93 are some constants. Its solution is N{t) = p{t) that is the Weierstrass elliptic function. 
The corresponding EoS is given by 

p=-p-12p\t)+g2. (9.2) 

The /-versions of Eq.(9.1) are /^ = if - g^f - 33, Jr = 4/^ - 52/ - 33 and /^ = - 32/ - 53- 
Another interesting models follow from the Bernoulli equation: 

N ^qi{t)N" +q2{t)N. (9.3) 

For example, if n = 2,gi = —t^,q2 — 2t^^ then the Bernoulli equation has the solution N — 
t'^{0.2t^ + const)-^. The /-versions of this solution are /(i?) = R^{Q.2R^ + const)-'^ and /(G) = 
G^(0.2G^ -I- const)~^. Finally we would like to note that it would be interesting to make relation 
with viable models of F{R) and F{G) gravities unifying inflation with dark energy (see e.g. [55] - 
[38] ) . This important question will be the subject of the separate investigation. 
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